Dynamic Jahn- Teller Mechanism of Superconductivity in MgB2 
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We propose a novel mechanism of superconductivity in MgB2 based on the dynamic electronic 
structure of the pcr-orbitals coupled with e2g phonons. A nonconventional superconducting state 
is found to arise from electron-phonon interactions in the presence of additional pairing channels 
made available by the dynamic Jahn- Teller effects. A partially broken pseudo-spin symmetry in 
this Jahn- Teller system, together with two-phonon exchange pairing, naturally gives rise to two 
distinct gaps both of which are basically isotropic in the {kx, ky) space. Important experimental 
observations including high and the anomalous specific heat are explained using this theory. 
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PACS numbers: 74.20.Fg, 71.70.Ej, 74.25.Jb, 74.70.Ad 

Following the discovery of high Tc superconductivity 
near 40 K in MgB2 [Q, a lot of efforts have been de- 
voted to revealing the underlying mechanism of this novel 
superconductor Now it is widely agreed upon that 
a relatively strong electron-phonon coupling is responsi- 
ble for the superconductivity of this material and there 
exit, at least phenomenologically, two superconducting 
gaps i, |, i |, |, 11 . It has been proposed 0, that 
the complex Fermi surfaces having both 2 dimensional 
(2D) cylindrical sheets and 3 dimensional (3D) tubular 
networks produce two gaps. However, it is still an open 
question whether one gap is mainly associated with the 
2D cylindrical Fermi surfaces and the other is associated 
with the 3D Fermi surfaces. 

In this letter, we propose a mechanism of supercon- 
ductivity in MgB2 based on a dynamic Jahn- Teller effect 
arising from the interplay between the doubly degener- 
ate pa electronic states and the e2g phonon modes. The 
hopping motion of holes of the pa character on the Boron 
layers is constrained by the 'pseudo-spin', which will be 
defined later, and partial breakdown of the pseudo-spin 
symmetry is also important in the pairing mechanism. 
Unlike existing theories, we will explain both high Tc 
and two distinct gaps in MgB2 without invoking the ad- 
ditional 3D Fermi surface of pir character. 

According to the band structure calculations ||l^, 
the pa states on the Boron layer form doubly degen- 
erate Eg bands at F (to be precise, along the FA line 
on which the component of k parallel with the Boron 
plane, k||, is zero), where we have small hole pock- 
ets presumably responsible for the superconductivity in 
MgB2. (Splitting of the band degeneracy away from F 
will be taken into account later.) These hole states are 
known to be strongly coupled with e2g phonon modes 
with fiuj « lOmeV |jl|, |l|. The Eg hole doublet coupled 
with e2g phonon doublet states on the Boron layer can 
be represented by the so-called E ®e Jahn- Teller model. 
A typical adiabatic potential energy surface of the E^e 
Jahn- Teller model is illustrated in Fig. 1, where qs and 



correspond to the two phonon coordinates. 

It is well known that one must take account of the 
dynamical phases of the holes and phonons in order to 
describe the E ®e system properly. When the ionic con- 
figuration traverses around the potential well in Fig. 1 
once, the corresponding orbital wavefunction of the hole 
(for simplicity of discussion, we leave aside the real spin 
wavefunction which simply follows the elementary BCS 
theory) acquires a phase of —1, as if the orbital part had 
spin 1/2. We call this a pseudo-spin of the hole and de- 
note its projected value by a = ±\. This subject was 
extensively discussed in the literature [|l6j in the context 
of the Jahn- Teller effect or Berry's phase. We also note 
that the dynamic Jahn- Teller effect in doped Cgo super- 
conductors was studied before, although the methods and 
results there were different from our work JTzt . 

We introduce a single site Jahn- Teller Hamiltonian per 
unit cell using mass-weighted coordinate for the sym- 
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FIG. 1: Adiabatic potential surface of the E ® e Jahn- Teller 
system and the doubly degenerate pseudo-spin states ( 1 1 ) , 
I - i>) of the hole. 
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FIG. 2: Diagrams for pairing interactions with phonon ex- 
change : (a) one-phonon exchange, and (b),(c) two-phonon 
exchange. The numbers indicate pseudo-spins. The obvious 
real spin indices are omitted for notational clarity. Diagrams 
with all the pseudo-spin signs switched are allowed as well 
(not shown). 

metrized displacement Qjij = 0,e), i.e. = y^M^Q^, 

ho^ ]^{vl + pI) + ]^^l{ql + ql) + A{q0(j.:, + q^ay), (1) 

where A is the hole-phonon coupling strength and Ux , (Jy 
are Pauli spin matrices for the pseudo-spin representa- 
tion of Eg hole doublet states. In an analogy to the 2D 
harmonic oscillator, the vibrational mode in the (ge, q^)- 
space belongs to the eigenmode of the "phonon" angular 
momentum operator Lz {=qgpe — qePe) with the eigen- 
value of I, and the principal quantum number v becomes 
v — I + 2k (fc is a non-negative integer and —ly < I < ly). 
In the present theory of superconductivity, only v = 0, 1, 
2 and / = 0,±1 are relevant. Defining the total pseudo- 
spin of the system by Jz ^ + \<Jz, it is straightfor- 
ward to prove [ho, Jz] = 0, that is, the total pseudo-spin 
of the system is conserved, while individual Lz and cr^ 
are not. (The phonon angular momentum component, 
Lz , can be regarded as a part of the total pseudo-spin in 
this context.) These pseudo-spin states are derived from 
the E ® e Jahn- Teller coupling between pa holes and 
e-2g phonons, and reflect the non-adiabatic nature of the 
coupled states. As a consequence, relevant electron and 
phonon states are expected to be strongly renormalized 
and contribute to anomalous behaviors (e.g., broad line- 
widths) observed in recent photoemission spectroscopy 
and Raman experiments H, ||, [ij] . 

We generalize the theory to the Bloch state in the k- 
space of a crystal by linear combinations of the local holes 
of the same pseudo-spin, but mathematical details will be 
reserved for a longer paper . We first consider the one- 
phonon [v ^ \) exchange for the Cooper pair formation. 
Each phonon has Z = ±1 in this case. In the one-phonon 
absorption or emission, the pseudo-spin of the hole should 
flip to conserve the total pseudo-spin. The only allowed 
one-phonon exchange mechanism is the opposite pseudo- 
spin Cooper pair as illustrated in Fig. 2(a). Pairing of the 
same pseudo-spin holes violates the conservation law in 
either one of the two (emission and absorption) processes. 
Such a selection rule has already been known in the self- 
energy correction to the Jahn- Teller system Q . 

We next consider the two-phonon [v — 2) case. 
In MgB2, the two-phonon exchange process has been 



thought to be important in producing high Tc and in- 
cluded in some previous works ]Tot . When we properly 
consider the pseudo-spin symmetry of the Jahn- Teller 
system such as MgB2, incorporating the two-phonon ex- 
change turns out to be essential in order to obtain two 
distinct gaps. Unlike the one-phonon exchange process, 
both opposite- and same- pseudo-spin pairs are allowed 
here as illustrated in Fig. 2(b) and (c). The total I (sum of 
two individual Vs) of intermediate phonons is restricted 
to zero (i.e., the hole pseudo-spin change of ±2 is not 
possible) to conserve the total pseudo-spin. A new model 
Hamiltonian with different pseudo-spin channels explic- 
itly labelled is, in the k-space representation, 

n^no + Hi + H2 + n3. (2) 

Ti-o = ^ekcJ^k^Cakcr, (3a) 

akcr 

Hi = - XI ^kk'clak'Tci-k'iC-"-kiCakT, (3b) 
akk' 

H2 = - X! ^kk'Cak'T^-a-k'i^-"-kiC„kT, (3c) 

akk' 

H3 = - X Vkk'clk'TcLk'iCo-kiCokT- (3d) 

akk' 

Hi, 7^2, and Ha correspond to the diagrams in Fig. 2(a), 
(b), and (c), respectively. The explicit form of the derived 
interactions, V^y^,'s {i —0, 1, 2), can be calculated and 
will be given elsewhere To discuss the nature of 

superconducting states, it suffices here to treat them as 
given parameters. The additional pairing channels due 
to the presence of pseudo-spin symmetry can give rise 
to a non-trivial superconducting ground state. Although 
the pairing is still based on electron-phonon interactions, 
the pairing order parameter can now attain more than 
one symmetry representation leading to the ground state 
with multi-dimensional order parameters while retaining 
the conventional s-wave type superconductivity. 

In the pairing Hamiltonian described above, however, 
important features of a real crystal are missing. First, 
the pseudo-spin symmetry is exact only at k|| ==0 and 
there occurs small but finite deviation from the exact 
symmetry for k|| 7^ where Fermi surfaces actually ex- 
ist. The band splitting at k|| 7^ is a clear manifestation 
of the breakdown of the rigorous Jahn- Teller theory for 
Bloch wavefunctions of k| | 7^ in a crystal as opposed to 
the case of a highly symmetric molecule. Second, higher- 
order Jahn- Teller coupling terms (not considered so far) 
do not commute with Jz and break the pseudo-spin con- 
servation II21]] . A simple way to incorporate the symmetry 
breaking effect is to introduce a small symmetry breaking 
parameter (/)k (0k||=o = and ^kn^^o 7^ 0) in the hopping 
Hamiltonian (Eq. (3a)) as follows. 

Ho = X ffkC^kaCaktr - X 0k£kc],k„C_c«ka- (4) 
akff akcr 
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This replacement naturally introduces the symmetry- 
breaking hopping process which gives the band splitting 
away from k|| =0. (If the pseudo-spin i were continu- 
ous symmetry {SU{2)) for all k's, there would be a sin- 
gle finite gap plus a zero-gap Goldstone mode, unlike 
the experimental observation in MgB2.) Now we choose 
the following two order parameters, one for the opposite 
pseudo-spin pair and the other for the same pseudo-spin 
pair. 



A 



(5a) 
(5b) 



It is to be noted that the pseudo-spin pairing assumed 
above is a triplet state (even for exchange of two parti- 
cles) to preserve the singlet (odd) state of the real spin 
pair and the s-wave (even) state of the orbital wavefunc- 
tion. Following the standard mean field approximation, 
we can derive the effective Hamiltonian H' (= H — epJ^), 



n' 



E 



^k(c^kT'^"kT C_Q_k|cLQ,_ki) 



- ^k(cj,ktC-akT - C-a-klcl_^^) 
ak 

- J2 ^l,k(cIkTcia-ki + C-a-kiC^kt) 

Qfk 

- A2,k(c|;,kTcLki + Ca-kiCakt) 
ak 

+ (constant terms), 



(6) 



where Ck = £k-£F, <^k 
ters are expressed as 



'kEk, and the gap- like pareme- 
Ai,k = 5](VkV + VkV)Ak', A2,k = 5]VkVi?k'. (7) 



Using the basis (ct j^^, ij^^, ci _kii c_ i _kx)i it is con- 
venient to cast the Hamiltonian into a matrix form (by 
omitting of the subscript k for notational simplicity), 



-S 

s e 

A2 -Ai 
V-Ai -A2 



-A2 -Ai\ 
-Ai -A2 

-e s 



(8) 



This Hamiltonian is easily diagonalized, for each k, yield- 
ing two positive eigenvalues. 



i?± = V(eT<5)2 + (A2±Ai)2. 



(9) 



The larger gap A+(= A2 + Ai) is associated with the 
band of i?normai = £, — S and the smaller gap A_ {— | A2 — 
All) is associated with the band of -Enormai — £, + S. 
Because of this band splitting, two gaps are associated 
with two different normal-state (above Tc) densities of 
states (DOSs). If we define N to be the normal-state 



DOS of unsplitted pa band at the Fermi level before the 
symmetry-breaking (j) is introduced, normal-state DOSs 
of two bands are N± — with the assumption of (f>k's 
being independent of k (^k = 0) in the vicinity of the 
Fermi level. Experimentally, it is observed that the state 
involved in the larger gap has the larger normal-state 
DOS than that of the smaller one j^]. In this situation, 
we naturally assign a positive value for (iV+ = > 

We first follow, for heuristic purposes, the approxima- 
tion procedure of the weak-coupling ECS theory. We 
assume some effective cutoff energy lo and define Ai = 
N{V° + V^) and A2 = NV"^. Only Ai, A2 and uj are in- 
dependent variables for the gap equations. The coupled 
gap equations for finite temperature read as 



A2 + Ai 

l-4> 
. A2 



Ai 



A_ = 



l + ( 
A2 + Ai 



A2 — Ai 
' l-ch 



— — tanh — — — 
2i?+ 2 

(j^ ^ „ tanh ^ 



•2E^ 

d£ tanh — 

' A+ 
d^^^tanh 
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(10) 



By fitting to Tc and the specific heat behavior, we are 
going to obtain below Ai = 0.235, A2 = 0.25 and (f> = 0.2 
for the assumed value of huj = lA.bmeV. The normalized 
solutions of the above gap equations (with ~ 0.94 and 
(j> = 0.2) are shown in Fig. 3(a). Our theory gives two dis- 
tinct gaps of the identical Tc originated from two split-off 
pa bands alone, in contrast to other theories which invoke 
the pir band to obtain multi-gaps. Furthermore, we pre- 
dict that there exist two gaps for the same k vector and 
both gaps are basically isotropic in the Boron plane (ex- 
cept for a small anisotropic influence from the pn bands 
located at different regions of the k space) unlike other 
theories proposing generically anisotropic gaps. With- 
out pseudo-spins, there would be only one s-wave gap 
solution giving the lowest free energy. It should be pos- 
sible, with angle-resolved experimental probes, to test 
the symmetry or angular dependence of the gap order 
parameters. Also, it is important in the present theory 
to include the symmetry-breaking parameter (/)k; other- 
wise, two gaps should be exactly degenerate and only one 
gap would manifest itself. With the same degree of ap- 
proximation as in the BCS formalism, we can derive the 
expression for Tc (weak coupling limit). 



ksTc ~ 1.13fi.wexp 



(Ai + A2) - V(^i-^2)2+4(/.2AiA2 



2A1A2 



(11) 



We then calculate the specific heat with this gap struc- 
ture. The resulting specific heat is plotted in a usual 



4 





1 - 


o 




II 




+ 

< 


0.5 - 


p 




< 


- 


1- 


1.5 - 




1 - 






CO 
CD 

o 


0.5 - 
- 





FIG. 3: (a) Solution of the self-consistent equations for two 
gaps A+ and A_, and (b) calculated specific heat for the 
above gap structure. 



manner in Fig. 3(b). To obtain both Tc = 39 K and 
the specific heat data in agreement with experiment, we 
require (f) — 0.2, Ai = 0.235, A2 = 0.25. The symmetry 
breaking parameter </> — 0.2 gives about the same band 
sphtting as ab initio calculations for the two pa bands 
near the Fermi level. Needless to say, since Tc is high 
and the material parameters are not in the weak-coupling 
limit, Eq. (11) is not to be applied directly to MgB2. We 
have to use a more realistic (Eliashberg) equation which 



1.20 '^-'^P 



1.04(1+A) 



gives, for example, Tc _ t^o" ^^i' [~ a-(i+o.62A)m» 
It is beyond the scope of the present paper to calculate 
Ai and A2 from the first principles. The effective cou- 
pling constant A in the exponent in Eq. (|l^), namely, 
j2Mi should be about 0.95, as is 

(Ai + A2)-V(^i -^2)^+40^ A1A2 

well known ||l3|, |l^, in order to give Tc — 39 K with 
H* = 0.13. Then wc obtain Ai = 0.71, A2 = 0.76 by sim- 
ple scaling of the previous values. These values are not to 
be taken literally because the two-phonon exchange can 
modify the empirical expression for Tc of the Eliashberg 
equation. In any cases, two-phonon exchange together 
with extra pairing channels makes an important contri- 
bution to Tc and it is the source of larger effective A in 
MgB2 than in other materials. 

In summary, we have proposed a pairing mechanism in 
superconducting MgB2 based on the dynamic Jahn- Teller 
effect of the doubly degenerate pa hole states coupled to 
degenerate 623 phonons. The pseudo-spin contraint im- 
posed on the motion of holes by the Jahn- Teller inter- 
actions is an important factor in understanding the dy- 
namical properties of the system. The one-phonon and 



two-phonon couplings lead to effective pairing interac- 
tions with additional pseudo-spin channels. As results, 
the high Tc observed in MgB2 is attributed mainly to the 
two-phonon exchange mechanism, and further the sym- 
metry breaking hopping term plays a crucial role in pro- 
ducing two gaps which are responsible for the anomalous 
behaviors observed in experiment including the specific 
heat measurement. 
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